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Abstract. We review recent developments in the field of microscopic transport 
model calculations for ultrarelativistic heavy ion collisions. In particular, we 
focus on the strangeness production, for example, the 0-meson and its role as 
a messenger of the early phase of the system evolution. Moreover, we discuss the 
important effects of the (soft) field properties on the multiparticle system. We 
outline some current problems of the models as well as possible solutions to them. 



1. Introduction 

The purpose of this overview is to present recent developments in the field of 
microscopic transport theory for ultrarelativistic heavy ion collisions, that is, for 
energies at SPS and RHIC. Special emphasis is put on the strangeness production. 
Moreover, we will particularly focus on the role of the (/>-meson as well as on the 
importance of the color field properties on strange observables. Furthermore, we will 
also briefly address the topics baryon transport, excitation functions, multi-particle 
collisions, initial conditions, J data, and particle interferometry. 

2. Color Fields 

An example of what is meant by the importance of field effects on the dynamics 
is given by the quark molecular dynamics (qMD) model |B G] ■ It is use d to study 
microscopically the dynamics of colored quarks. The intention is to get an effective 
description of the non-perturbative, soft gluonic part of QCD. In particular, the 
hadronization process itself can be described explicitly. The initial quark distributions 
are obtained, for example, from the UrQMD model. The non-equilibrium dynamics of 
hadronization and the loss of correlations among quarks can be studied subsequently 
by qMD. The semi-classical model has a two-body color potential between the 
(anti)quarks. In addition, one has to define a dynamical hadronization criterion. The 
quarks interact via a Cornell potential with color matrices. Thus, they carry color, 
(flavor, spin, and isospin). The simple model Hamiltonian is 

JV 

H 

*=1 ~ i,j 

§ ssoff@th.physik.uni-frankfurt.de 



N 



News on Strangeness at Ultrarelativistic Energies - Review of Microscopic Models 2 



where N are the number of quarks in the system. The potential V(r) is linearly 
increasing at large distances (confinement) and has a Coulomb-type behavior at short 
distances. The Cornell-potential is 

. . 3 a s 

V [rj = \- nr . 

4 r 

A further approximation is done by allowing only the diagonal colorless gluon 
exchange, that is, the quarks always keep their color. The color matrix elements 
Cij determine the sign and relative strength of the interaction depending on the color 
combination of the pair (for details see P). 

The rather simple model is surprisingly able to describe a large amount of 
experimental data, for example, various particle spectra and ratios. Concerning 
strangeness, we show here the interesting feature of aquiring net strangeness in the 
quark phase. According to the idea of the strangeness distillation antistrange quarks 
hadronize earlier in a baryon rich system. Fig. 1 shows the number of net strange 
quarks (s-s) as a function of time for the baryon rich system in central Pb+Pb collisions 
at 30AGeV. Values up to approximately 10 are reached depending on the string 
tension. This indicates the possibility to study hyperon matter, hypernuclei, or even 
strangelets in these high baryon density collisions around the future GSI energies. 

Strangeness Distillation at Pb + Pb @ 30 GeV/N 




Figure 1: qMD calculation for 
Pb+Pb at 30ylGeV by S. Scherer 
pQ. The number of net strange 
quarks is shown as a function of 
time for different values of the 
string tension k. 
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Another feature of the model is to study the details of the parton-hadron phase 
transition that is mainly driven by quark rearrangement. Late stage quark-antiquark 
pair production is a rather rare process in this model because the strong color fields 
that are needed are screened. Fig. 2 shows the distribution of the mean path length of 
quarks before hadronizing. Quarks from the same original hadron (solid line) recluster 
faster than quarks from different original hadrons (dotted line) . The hadron formation 
follows an exponential decay of the quark cluster to hadrons. The reclustering exhibits 
a different decay length than the rearrangement process (2.2 vs. 4.8 fm). The amount 
of reclustering (ss 50% for S+Au) demonstrates the survival of correlations in the 
quark system and possibly incomplete thermalization. This kind of detailed analysis 
also provides additional insight to the question of the mechanisms of hadronization 
(recombination vs. fragmentation 0). 
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Figure 2: Quark Molecular Dy- 
namics simulation of the hadroniza- 
tion in S(200AGeV)Au collisions. 
Quarks originating from the same 
hadron hadronize after a shorter 
path length s than quarks stem- 
ing from different hadrons. Figure 
taken from |J. 
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3. </>-mesons, baryon transport and the sensitivity to field properties 

The field properties also enter in the description of heavy ion collisions by means 
of cascade calculations that are based on hadron, (di)quark, and string degrees of 
freedom, as for example the UrQMD model Here the initial particle production is 
dominated by string excitations and fragmentations, that is, longitudinally excited 
color flux tubes. In a dense colored medium the color field strength might be 
considerably enhanced. The particle formation process in hadronic collisions can be 
viewed as quantum tunneling of quark- antiquark and gluon pairs in the presence 
of a background color electric field. It is formed between two receding hadrons 
which are color charged by the exchange of soft gluons while colliding. In nucleus- 
nucleus collisions the color charges may be considerably greater than in nucleon- 
nucleon collisions due to the almost simultaneous interaction of several participating 
nucleons. With increasing energy of the target and projectile the number and 
density of strings grows, so that they start overlapping, forming clusters, which act 
as new effective sources for particle production |SJ EJ. It has been predicted that 
the multiplicities of, for example, strange baryons or antibaryons should be strongly 
enhanced |S| once the color field strength grows. The abundances of (multiply) 
strange (anti)baryons in central Pb+Pb collisions at Cern-SPS J2J, for example, can 
only be explained within the framework of microscopic model calculations ja] if the 
elementary production probability of ss pairs, which is governed in the string models 
|H] by the Schwinger mechanism ^U] ~ exp(— / nm 2 q /2n) 1 is considerably enhanced. 
This corresponds either to a dramatic enhancement of the string tension k (from the 
default ~ 1 GeV/fm to 3 GeV/fm) or to quark masses m q that are reduced from their 
constituent quark values to current quark values as motivated by chiral symmetry 
restoration. In a fully analog way, p abundances can be explained since an enhanced 
string tension similarly leads to an increased production probability of antidiquark- 
diquark pairs which is needed to account for the experimentally observed yields 
However, this argument, providing additional motivation, has to be reconsidered if 
multifusion processes e -g-> 5 71 " ~~ * BB, that are neglected in microscopic transport 
models based on 2 — ► n scatterings, contribute significantly to the baryon pair 
production. A variation of the string tension from k = 1 GeV/fm to 3 GeV/fm 
increases the pair production probability of strange quarks (compared to light quarks) 
from 7 S = P(s~s)/ P(qq) = 0.37 to 0.72. Similarly, the diquark production probability 
is enhanced from j qq — P(qqqq) / P{qq) — 0.093 to 0.45. In general, heavier flavors or 
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diquarks (Q) are suppressed according to the Schwinger formula ^U] by 



1Q 



P(QQ) 

P(qq) 



exp 



ir{m 2 Q - m|) 
2k 



The effectively enhanced string tension in a densely colored environment also 
follows from its relation to the Regge slope a'. In the rotational string picture the 
string tension k is related to the Regge slope a' by |131 HI) 

1 



2-ko! 



The empirical value of the Regge slope for baryons is ol w 1 GeV~ 2 ^5] that 
yields a string tension of approximately 1 GeV/fm. However, high-energetic processes 
dominated by Pomeron exchange, characterizing the multi-gluon exchange processes 
existent in high-energetic nucleus nucleus collisions, are described by a Regge 
trajectory (Pomeron) with a smaller slope of a' P « 0.4 GeV -2 ^JEl- According 
to Eq. @ this translates into a considerably larger {effective) string tension k. 

Strong color fields also have an impact on the formation times of newly produced 
(di)quarks and thus on collision rates. The formation times are inversely proportional 
to the string tension 

t f ~ 1/k. 

The larger the string tension the shorter and short-living are the strings with a certain 
total energy. 




(j) - rapidity 

Figure 3: Rapidity distributions of ^-mesons in central Au+Au collisions at RHIC. 
Besides the distributions of all tfi'a (full circles), the contributions from 0's produced 
in (resonant) meson-meson (KK) collisions (open circles) and from cj>'s originating 
from string decays (open squares) are shown. Each panel corresponds to different 
values of the string tension (k) (mass term) and the intrinsic transverse momentum 
broadening parameter a. 

Moreover, Schwinger's formula does not only apply to the mass term but to the 
transverse momentum term, too. While k determines the production probability, 
the intrinsic transverse momentum parameter a defines the transverse momentum 
scale. The vacuum value of 0.55 GeV/c (width of a Gaussian) may be considerably 
enhanced by the strong color fields. Here, we vary both parameters independently 
to study systematically their individual effects on the observables JH|- Fig- 3 shows 
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the rapidity spectra of </>-mesons at RHIC. The intrinsic p t is varied from 0.55 to 1.5 
GeV/c (left to right). The upper and lower panels correspond to the vacuum string 
tension and the SCF scenario, respectively. Clearly a larger k value leads to a strongly 
enhanced <p production. Increasing the intrinsic pt decreases the </>-meson yield. The 
reason for this becomes immediatley clear when looking at the production channels. 
The production in meson-meson (KK) collisions gets considerably reduced when tr 
is enhanced. The coalescence-like production loses due to a growing phase space. 
The experimental study of the </>-meson in both decay channels and the observed 
differences (due to rescattering effects) at RHIC will help to disentangle the different 
source contributions and thus enlight the role of strong color fields. 

Another observable that has been shown to be highly sensitive to SCFs is the 
netbaryon number transport. Prior to the collision it is non- vanishing only at beam 
and target rapidities (Ay ~ 11 units for RHIC!). Recent data by the STAR |19j . 
PHOBOS H3, BRAHMS [2T], and PHENIX |23 collaborations, show a considerable 
number of netprotons (p — p) at midrapidity. and antibaryon-to-baryon ratios smaller 
than unity. 

The calculated anti-baryon-to-baryon ratios at midrapidity are shown in Fig. 4 
as a function of the strangeness content both for k = 1 GeV/fm as well as for 
the strong color field scenario (SCF, k = 3GeV/fm) [23, ■ The £?/_B-ratios increase 
with the strangeness content \S\ of the baryons and also increase with with impact 
parameter b. There is stronger absorption of antibaryons in central collisions (reducing 
the numerator of the ratio) and more stopping (increasing the denominator). Most 
important, the i3/i3-ratios also increase with the color field strength k what is basically 
due to an increased pair production (while more collisions should reduce the ratios). 
The SCF results are closer to data. Also the observed hyperonization in the netbaryon 
distributions supports the SCF scenario (the (A — A)/(p — p) ratio is ~ 0.8 for the 
strong color field (SCF) (k = 3GeV/fm) scenario, corresponding to the experimental 
value [221 EH] but smaller for the vacuum calculations. 

Figure 4: Antibaryon-to-baryon ra- 
tios at midrapidity as a function 
of the strangeness content \S\ in 
Au+Au collisions at RHIC (i/s NN = 
200 GeV). Calculations with a string 
tension of n — lGeV/fm are shown 
on the left and the results with strong 
color fields (k = 3GeV/fm) are 
shown on the right [23] . 

A complementary or alternative ansatz to describe the baryon transport is 
provided through the parton cascade model (PCM). Fig. 5 shows the netbaryon 
distributions as calculated with the VNI/BMS version |^][23] of the parton cascade. 
It is interesting to note that a considerable fraction of the observed netbaryon number 
at midrapidity is already contained in the parton distribution functions of the original 
nuclei. Other contributions arise from secondary scatterings of the partons and 
fragmentation that together yield a netbaryon number compatible with data. It will be 
interesting to see whether the correct longitudinal dynamics will also imply the strong 
transverse expansion dynamics as, for example, the relatively large elliptic flow. 
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Figure 5: Netbaryon rapidity 
distributions for Au+Au at RHIC. 
Crosses are calculation with 
only primary parton scatterings, 
rhombes include parton rescatter- 
ing and squares include rescattering 
and parton fragmentation. The 
solid lines show the netbaryon 
content of the partonic distribution 
functions for gold nuclei, scaled by 
an average liberation factor of 0.4. 
Figure taken from |24| . 
The production of 4> mesons in Au+Au collisions at RHIC has also been studied 
with the quark gluon string model (QGSM) 26 . It is found that the inverse slope 
parameter agree with data but the absolute yield of </>'s is underestimated by a factor 
2. Somehow surprising, the fusion of strings does not increase the (j) yield in QGSM. 
Similar to previous UrQMD results 8 the early production is governed by string 
processes. The resonant kaon rescattering production channel dominates in the later 
stages. Fig. 6 illustrates the freeze-out distribution of the mesons over rapidity and 
emission time from QGSM. 
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Figure 6: d 2 N/dydt distribution 
of the final-state <f> mesons in the 
(y, t)-plane from QGSM ^ 



The puzzle has also been addressed with a systematic study of the kaon and 
dimuon decay channels by the AMPT model [57|. The AMPT model is a multiphase 
transport model that includes the initial partonic and final hadronic interactions. It 
starts with minijet partons and strings from the HIJING mode. The partons enter the 
ZPC parton cascade model. The parton hadron transition is based on the Lund string 
fragmentation model. Final-state hadronic scatterings are modeled by the ART model. 
Similar to the other models the resonant KK production channel also contributes 
considerably to the total <j> yield. However, even with in-medium modifications of the 
kaons and <fi's the observed differences in the yields cannot be fully explained (see 
figure 7). An additional ad hoc production is needed to come close to the data (of 
NA50) . It is important to remark that the initial (j>'s (those that are possibly created 
through strong color fields) will be primarily visible in the dimuon channel. 
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Figure 7: Transverse mass spectra 
of (j>'s reconstructed from K + K~ 
pairs (solid symbols) and from n + u~ 
channel (open symbols) for Pb+Pb 
at 158A GeV {b < 3.5, \y\ < 
1) fm in the AMPT model. The 
results are for without (squares) 
and with (triangles) in-medium mass 
modifications. The solid circles are 
NA49 data gH] for <j> -> K+ K~ 
decay, and the open circles are NA50 
data |29| for <f> — ■> In the 

inset is shown as a function of raj- the 
ratio R(niT) for phi mesons decaying 
to kaon-antikaon pairs that are not 
scattered to those determined from 

1 8 

' the dimuon channel. Circles are for 
an artificially increased number of <^>'s 
by a factor 2 in the HIJING model. 
Figure from |27|. 



4. Further Topics 

Excitation functions of meson abundancies and multi-meson fusion 

New data from the Cern/SPS allow us to study the excitation function of particle 
production (and spectra, see (30)) m the interesting energy range of maximum 
netbaryon densities. Characteristic features in the data as maxima in the K + /n + 
ratio gave reason to speculate about the observation of the phase transition. How 
do 'conventional' hadron transport theories compare to data [31 1 1321 I33|? The 
detailed comparison with two realisations of the cascade model with hadron resonance, 
(di)quark, and string degrees of freedom (UrQMD and HSD) shows that the kaon 
yields are described reasonably (see Fig. 8). However, the pion yields and hence the 
K+/-7T+ ratios are off (see figure 9). 
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Figure 8: Excitation function of K + and K yields from central Au+Au (AGS) or 
Pb+Pb (SPS) collisions (HI) in comparison to experimental data 135) . 

The overestimation of the pion yields in the model calculations could be due 
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to a lack of multiparticle back reactions in the model that also lead to 'shorter' 
chemical equilibration times. 2 — > n collisions (with n > 2) are performed but the 
inverse channels are not taken care of. Another possible explanation is given through 
enhanced pion in-medium masses. In fact, the importance of 3 <-> 2 transitions has 
been explored in an extended HSD transport approach for antiproton production 
by meson fusion in A + A collisions at the AGS and SPS (see figure 10). In order 
to achieve a more conclusive answer from transport studies multi-particle interactions 
deserve further investigation in future generations of transport codes. 
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Figure 9: The same as in Fig. 8 for the K + /tt + and K jn ratios. 




Figure 10: Illustration of the flavor rearrangement model for BB annihilation to 3 
mesons and vice versa. The mesons Mi may be either pseudo-scalar or vector mesons, 
respectively. Figure from |36) 

Within this approach the abundancies of antiprotons as observed from peripheral 
to central collisions of Pb + Pb at the SPS and Au + Au at the AGS can approximately 
be described. 

Fig. 11 shows the reaction rate B + B — > mesons vs. the backward reaction rate 
for central Pb+Pb collisions at 160 AGeV. Both rates turn out to be comparable 
within the statistics demonstrating an approximate local chemical equilibrium. The 
final p's origin dominantly (~ 88%) from the 3 meson fusion reactions. These studies 
indicate that these back reactions may have dramatic effects and therefore ask for 
further investigations. 
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Figure 11: The annihilation rate 
BB — > mesons (dashed histogram) 
for central Pb + Pb collisions at 
160 A GeV as a function of time in 
comparison to the backward reac- 
tion rate (solid histogram) within 
the HSD transport approach. Fig- 
ure from |36| . 
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4.2. J/tp 

Another possible QGP signal is charmonium production and its anomalous 
suppression. It is worth to mention that the new data from NA50 however, do 
not deviate significantly from the predictions of standard-type cascade calculations 
(UrQMD or HSD) jS3 EHj (Fig. 12). Both the E T dependence and the E ZDC 
dependence do no longer support the observation of a strong anomalous suppression. 




Figure 12: The ratio of J/\P to Drell-Yan production as a function of forward 
energy Eznc (left) and transverse energy Et (right) for Pb+Pb at 158 AGeV. The 
preliminary data from the year 2000 NA50 Pb+Pb run |.S9) are compared to UrQMD 
results |37) . Figure from |3.S) 



4-3. Initial conditions 

The newly developed NEXUS3 model addresses the initial conditions of pp or AA 
collisions 001^2- It provides a self-consistent description of the energy sharing among 
the participating partons in the collision. The idea is called parton based Gribov-Regge 
theory, that is, soft and hard processes are described by (multi-)pomeron exchange. 
As one consequence strings are not necessarily attached to the valence quarks but may 
be connected to sea quarks. An observable supporting this picture may be given by 
the fj /CI ratio in pp collisions which is smaller or close to unity for NEXUS3 but larger 
for standard string models (or the older versions of NEXUS) |42j . 
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4-4- Particle inter •ferometry 

In brief, progress has been made at several frontiers. More realistic models were 
developed and used to calculate the multidimensional correlation functions. Hybrid- 
type models, for example, that allow one to study explicitly the sensitivities to the 
QCD equation of state by simultaneously describing the freeze-out in a realistic way 
were used to discuss the so-called RHIC HBT-puzzle. The impact of in-medium 
modifications on HBT results were discussed for the first time. For a discussion and 
references, see for example (431 1441 EES] . 

5. Outlook 

Topics not discussed here are, for example, the necessary tests of the microscopic 
models by simulating infinite matter. The upcoming RHIC pA data will put severe 
constraints for the models. An exact understanding of the elliptic flow V2 is needed 
(collision numbers, formation times). How is the (partially) hydrodynamic behavior 
realized on a microscopic level? Also, there has been quite some progress in the charm 
sector, for example, the prediction of a substantial production of J/v&'s in DD collisions 
|4"fi] , in analogy to the 4> production via KK [H] . 

6. Summary 

Here, we focused on the role of strong color fields that lead to (i) enhanced heavy 
flavor and diquark production probabilities, (ii) modified, shorter formation times, 
(iii) stronger intrinsic transverse momenta. The effects on the <fi meson include (a) 
enhanced yields, (b) a change of the composition of the spectra, (c) a hardening of 
the spectra, (d) modified baryon dynamics (hyperonization). Moreover, we briefly 
discussed the competition of late stage vs. early production and the rescattering of 
the decay daughters. The differences between the kaonic and muonic ^-mesons are 
desired! 

The quark molecular dynamics model provides a (simple) dynamical description 
of the hadronization process in the soft regime and needs to be further tested. 

Parton cascades are now technically established and provide a valuable tool to 
learn about the validity range of pQCD or necessary modifications. 

The study of excitation functions was shown to be essential to detect deficiencies 
in the models. Additional new developments include the investigation of multi-particle 
collisions, of charm in microscopic models, and the formulation of the initial conditions 
in a field theoretical language (parton ladders) . 
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